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ABSTRACT 

We have used the AAOMEGA spectrograph to obtain R ~ 1500 spectra of 714 stars that are members of 
two red clumps in the Plaut Window Galactic bulge field {I, b) = 0°, —8°. We discern no difference between 
the clump populations based on radial velocities or abundances measured from the Mgb index. The velocity 
dispersion has a strong trend with Mgfe-index metallicity, in the sense of a declining velocity dispersion at 
higher metallicity. We also find a strong trend in mean radial velocity with abundance. Our red clump sample 
shows distinctly different kinematics for stars with [Fe/H] < — 1, which may plausibly be attributable to a 
minority classical bulge or inner halo population. The transition between the two groups is smooth. The 
chemo-dynamical properties of our sample are reminiscent of those of the Milky Way globular cluster system. 
If correct, this argues for no bulge/halo dichotomy and a relatively rapid star formation history. Large surveys 
of the composition and kinematics of the bulge clump and red giant branch are needed to define further these 
trends. 

Subject headings: Galaxy: formation — Galaxy: bulge — Galaxy: kinematics and dynamics 



1. INTRODUCTION 



The existence of a prominent bar in the Galactic bulge is 
now well e stablished from multipl e lines of e v idence (e.g., 
Liszt & Burton 1980; BHtz & Spergil [1991t IStanek et all 



1994t Babusiaux & Gilmore 1995). The large-scale kinemat- 



ics of the bulge sa mpled b y the BRAVA (Bulge Radial Veloc- 
ity Assay) survey (iRich et al. 2007; Howard et al. 2008, 2009) 
can be fitted with simple cylindrical rotation and little or no 
classical spheroidal component (Shen et al. 2010), a conclu- 
sion also reac hed (albeit at lower con fidence) from proper mo - 
tion surveys (iRattenbury et al.ir2007D . Our 'bulge' appears to 
consist of a pean ut-shaped bar, reminiscent of t he 'pseudob- 
ulges' defined by Kormendv & Kennicutt ('2004') and encoun- 
tered elsewhere (Kormendv & Barentine 2010). 

At the same time, this result is puzzling. The BRAVA 
data imply that our Galaxy has not undergone any significant 
merger since the epoch at which the disk formed, in contrast 
with expectations from simulations within A C DM cosmo- 
logical models (e.g.. Bullock & Johnston 2005; C ooper et al.l 
|2010). However, it is well known that stars in the bulge are 
metal-rich and a-enhanced, indicating a rapid sta r formation 
history typical of c lassical bulges (Mc William & Rich 1994; 
iBallero et al.ll2007h. The observa tion of an abundance gradi- 
ent for bulge stars (IZoccali et al.l l2008) and a correlatio n be- 
tween abundance and kinematics (Babusiaux et al."2010) may 
also indicate the presence of a classical spheroidal compo- 
nent, although dynamical data dShen et al.ll2010l) seem to be 
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inconsistent with this . On the other hand, detailed abun- 
dances for subgiants and lensed dwarfs, are more similar to 
those of the thick disk, in accordance with a pseud obulge for- 
mation scenario ( Melendez et al. 2008; Alves-Brit o. A. et al] 
IMoHBensbveTaLIGOlftlRvde et al.li201 0). 

Adding to the complexity, the red clump in the 

bulge col or-magnitude diagrarn appears doub led at 

\l\ > 5° j Mc Willi am & Zoccalil lIoTft iNatafetalJ IMoh . 
iMc William & Zoccali! (|201^ carry out a very careful analysis 
of the photometric properties of the double red clump and 
conclude that the splitting of the red clump is due to a 
distance effect and that the bulge may contain an X-shaped 
structure, extending from the ends of the bar So far, the 
only sp ectroscopic study of thi s population has been carried 
out by iRangwala et al.l (120091) using Fabry-Perot imaging 
spectroscopy in Baade's Window and some adjacent fields, 
and their results suggest that there is a metallicity gradient 
with galactic latitude (Rangwala & Williams 2009), as well 
as dynamical differences, which would be more consistent 
with the presence of a stream, although it is clear that more 
accurate kinematics and metallicities are needed to truly 
understand the nature of this feature. 

Here we present the first low resolution survey of a field at 
I = 0° b = —8° (Plaut's Window) specifically targeting the 
red clump stars. We analyze the kinematics and metallicity 
of stars belonging to each peak and derive the metallicity- 
kinematics trends. Although this is clearly a 'first look' ex- 
ploratory analysis, our data imply that the two populations are 
dynamically and chemically similar and fav or the X -shaped 
bulge hypothesis proposed by Mc William & Zoccalil (I2bl0i) . 

2. DATA REDUCTION AND ANALYSIS 

The data for this project were kindly prov ided by the 
Galaxy and Mass Assembly Survey (GAMA - iDriver et al.l 
I20T0I) . This survey is observing three 12 deg^ fields at 9h, 
12h and 15h RA using the AAfl multi-fiber spectrograph on 
the Anglo- Australian Telescope (AAT). At the end of May, 
2010, it was found to be impossible to reach either GAMA 
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Fig. 1 . — The color-magnitude diagram (from 2MASS data) for our bulge 
field is shown in the left-hand panel of this figure. The blue box shows the 
(non-dereddened) selection limits used. The right-hand panel shows the his- 
togram of the distribution of stars in Kq and the two peaks identified by the 
GMM algorithm, marked by arrows (red for the blighter peak and green for 
the fainter peak). 

field during the last two hours of the night, at sufficiently low 
airmass. The survey kindly offered to observe one of the fields 
containing a double red clump for us. 

Two 400-fiber configurations were observed for Plaut's 
Window (where the two red clumps are relatively well sep- 
arated and the extinction is lower ). The target st ars were se- 
lected from the 2MASS survey (ISkrutskie et al.l 12006), in a 
2 degree (diameter) field. Bulge stars were required to have 
Kq = 7.5{J — K)o + 9to exclude disk contamination and to 
lie between 12.5 < Kq < 13.5 and 0.58 < ( J - K)o < 1.00 
to probe th e double red clump. Magnitudes were dereddened 
from the Schle gel et all (fT998! ) maps. Figure 1 shows the 
color-magnitude diagram of the target stars. The two red 
clumps are identified by appl ying a Gaussian Mixtur e Mod- 
elling algorithm to the data dMuratov & Gnedinll2010h , which 
also assigns to each star a probability of belonging to either 
peak. The magnitudes of the two red clumps peaks are in 
good agreement with those reported by McWilliam & Z occalil 
(|2010) for this region. By necessity, we used the same spec- 
troscopic setup as employed for the GAMA survey; this cov- 
ers the entire optical window, from about 3700 to 9000 A 
at a resolution of about 1500. Although this is somewhat 
less than optimal for determining stellar radial velocities and 
abundances, it suffices for our initial analysis of this field. 

The data were reduced using the automated pipeline sup- 
plied by AAil. A total of 714 stars were eventually observed. 
The typical signal to noise ratios of these spectra are about 10 
at Ca H & K, 20 in the Mg6 region and 50 for the Calcium 
Triplet wavelength range. Radial velocities were measured 
by cross-correlation against stellar templates using the runz 
program (Saunders et al. 2004) which is specifically written 
for analysis of AAil data. Among the several templates avail- 
able in runz we imposed the choice of a K-giant template, 
if this was not automatically selected by the program. Of the 
714 targets, 63 1 returned a valid radial velocity ( > 95% prob- 



ability that the measured velocity is correct from the height 
of the correlation peak). Typical velocity errors are 1/10 of a 
resolution element, or 25 km s^^ for the R ~ 1500 of GAMA 
data, but no stars with a velocity error above 2cr (50 km s~^) 
were used for our analysis. 

We experimented with a number of techniques to measure 
metal abundances, including the non-SEGUE stellar parame- 
ters pipeline (e.g., iLi et al.|[2010l and references therein) and 
the calcium triplet (although red clump giants were outside 
the luminosity range of the calibrations). We eventually found 
that the most reliable measurements (when compared with 
the high res olution abundances measured for giants in Plaut's 
window by Johnson et al .ll20lT|) were given b v using the Mg& 
Lick index (Worthev 199 4l: IIbata & Gilrnorei ri995). We mea- 
sured this index and its errors (based on the CCD noise pa- 
rameters and the radial velocity errors) using the LECTOR 
softwar^l The typical error in measuring the Mg& index is 
about 5%. 

In order to derive metal abundances f rom the Mg& iiidex, w e 
used fitting functions for cool stars bv 'Wort hev et al.l (119941) . 
We assumed log g — 2.25 which is appropriat e for red clump 
giants based on the Padova isochrones (Marigo e t al.ll2008l) . 
Temperatures were estimated from the dereddened J — K 
color, using the T^f / - color rela tions by iH oudashel t et aP 
(120001) . This approach was used by Cote et al.i(.1999,) . among 
others, to measure the abundance of red giants in M31, al- 
though we found that only Mg6 can be reliably measured in 
our data. We derived metallicities for 545 stars in our data. 

The typical random error in metal abundance is ±0.1 dex, 
based on the error in the index measurement, but of course 
there are systematic errors depending on the assumed values 
of log .9 and T^ff. Altering \ogg by 0.1 at the same tem- 
perature yields a change in [Fe/H] of around 0.15 dex, while 
altering Teff by lOOK yields abundance changes of 0.3 dex. 
In addition, our data are not on the Lick system as no stan- 
dards were observed, and this may introduce a systematic ef- 
fect of the order of 0.2 dex in the measurement of metal abun- 
dances. We used som e Lick standards observed by ELODIE 
iMoultaka et alj (120041) to measure Mg6 and find that our abun- 
dances tend to be ~ 0.2 dex too low, which would bring our 
da ta in better agreemen t with high resolution measurements 
by iJohnson et al.l (12011). In addition, the known a enhance- 
ment of the bulge ("e.g- lMcWilham & Richl[T994l) means that 
the Mg6 abundance may overestimate the actual [Fe/H] of the 
stars, although it is probably a better proxy of the total metal 
abundance. 

3. THE NATURE OF THE RED CLUMP 

Figure 2 shows the radial velocity distribution for red clump 
stars derived from our data. The mean heliocentric radial ve- 
locity for the entire sample is —15 ± 5 km s^^ (rms) , which 
is equivalent to a velocity of —4 km s^^ (Galactocentric stan- 
dard of rest), consistent with data from the BRAVA survey. 
The radial velocity dispersion is 110 ± 5 km s^^, which is 
somewhat larger than what is measured for this galactic lati- 
tude by BRAVA, although we are looking at a fainter sample. 
Th e intrinsic dispersion an d its error were computed follow - 
ing Spaenh auer et aTl (1199 2). The higher velocity dispersion 
seems to be due to a low metallicity population that may not 
be present in the BRAVA data. 

Stars in the first (brighter) peak have a mean heliocen- 
tric velocity of —18 ± 6 km s^^ and velocity dispersion of 

http://www.iac.es/galeria/vazdekis/vazdekis_software.html 
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Fig. 2. — Radial Velocity distribution of stars belonging to each red clump 
peak (red for the brighter peak and green for the fainter peak as in Figure 1). 
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Fig. 3. — [Fe/H] distiibution of stars belonging to each red clump peak (red 
for the brighter peak and green for the fainter peak as in Figure 2). The black 
dashed histogram shows the distribution of metal abundances for red giants 
studied by Johnson et al. 20 11. 

109 ± 5 km s^^, while stars belonging to the second (fainter) 
red clump have mean velocity of — 6 ± 8 km s~^ and veloc- 
ity dispersion of 113 ± 9 km s^^. Within the errors, stars 
in both red clumps have the same mean velocity and veloc- 
ity dispersion and therefore appear to have the same kinemat- 
ics. This would appear to be at odds with claims for a differ- 
ence in the kinematics of the two red clumps (Ran gwala et al.l 
I2009h . although more sightlines are needed. Application of 
the Gaussian Mixture Modelling algorithm shows that there 
is < 50% chance that the velocity distribution is bimodal. A 
Kolmogorov-Smirnov two-sample test yields an 88% chance 
that the two distributions come from the same parent. The 
t-test and F-test show that the two distributions do not have 
significantly different means or variances. 

Figure 3 shows the distribution in metal abundance for 
stars in both clumps. The distribution for all stars has a 
mean [Fe/H] of —0.61 ± 0.03. The distribution is consis- 
ten t with that observe d for r ed giants in Baade's Window 
by iMc William & RichI (119941) . the high resolution metallic- 
ities in Plaut's window (overplotted in the figure) measured 
by Johnson et al. (201 1) and the expectations from the metal- 
licity gradient measured by Zoccali^t al. (2008), but with a 
systematic offset 0.2 dex) probably due to the lack of an 
absolute calibration based on Lick standards. 

Stars in the first peak have mean [Fe/H] of —0.55 ± 0.03 
with a dispersion of 0.58, while stars in the second peak have 



mean [Fe/H] of — 0.67±0.03 with a dispersion of 0.62. Again, 
the Gaussian mixture modeling algorithm returns a distribu- 
tion consistent with a unimodal distribution. Within the er- 
rors, stars belonging to each red clump have the same [Fe/H] 
abundances. A K-S test gives a 73% probability for the two 
clumps to be drawn from the same distribution. Similarly, the 
t-test and F-test show that two samples do not have signifi- 
cantly different means or variances. 

One aim in this work is to compare the kinematics and 
metallicities of stars belonging to the two red clumps iden- 
tified jnjhe galactic bulge by McWilliam &Zoccali (2010); 
iNataf et al.] (12010 ). Our data show that the two populations 
have the same kinematics and metal abundance, which sug- 
gests that the difference in luminosity between the two red 
clumps is a distance effect. In other words, the two red clumps 
are observed at the two ends of the ~ 2 kpc bar that consti- 
tutes the bulge of the Milky Way. These data are therefore 
consistent with the analysis by McWilliam &Zoccali (2010), 
attributing the double red clump to the existence of X-shaped 
protrusions at the end of the Milky Way bar (i.e., an X-shaped 
bulge). If these stars lie at the end of a bar, then the metal 
abundance gradient across the structure is expected to be quite 
small. 

An alternate possibility is that the two red clumps are dif- 
ferent in helium content. Helium abundance variations are 
believed to exist (from indir ect evidence) in ma ssive globular 
clusters in our Galaxy (e.g..'Carrett a et all200 9). Helium en- 
hancement would produce differences in the luminosity of the 
red clump stars (D'Antona et al . 2010), but any process capa- 
ble of producing the necessary helium enhancement would 
also overproduce metals. We would presume that populations 
with different chemical abundances might also have different 
kinematics. 

4. ABUNDANCE TRENDS WITH KINEMATICS 

Classical bulges present a number of abundance trends with 
kinematics , some of wh i ch are also observed in the Galac- 
tic bulge. IZoccah et al.l (|2008[) have measur ed a metallicity 
gradient with radius. iBabusiaux et al] (120101) find a trend be- 
tween metal abundance and velocity dispersion in Baade's 
Window and two lower (b — —6° and —8°) galactic latitude 
fields, arguing for a two component bulge, with a metal-rich 
system comprising the bar and a metal-poor spheroid or thick 
disk, the relative contribution from each of these varying with 
galactic latitude. 

We plot the mean velocity (lower panel) and velocity dis- 
persion (upper panel) as a function of metallicity in Figure 
4. We plot all stars (black points), stars in the first (brighter) 
peak (red points) and stars in the second (fainter) peak (green 
points) separately. Stars in each group are binned in bins con- 
taining the same number of stars after sorting by metal abun- 
dance. For each of the groups we consider (all, stars in the 
brighter and fainter peaks) the bins span a non-overlapping 
range in metal abundance. Stars in both red clumps appear 
to obey the same relations between metallicity and kinemat- 
ics. This confirms that the two populations do not differ in 
kinematic properties or abundances, which is more consistent 
with a projection effect and an X-shaped bulge structure. 

The data shown in Figure 4 show a clear trend for increasing 
velocity dispersion with decre asing metal abundanc e. This 
is similar to what is f ound by iJohnson et al.| ( 1201 ll) in their 
higher resolution data. IBabusia ux et al. (2010) uses two fields 
at 6 = —6° and b — —12° and although there are few stars 
at [Fe/H] < — 1 in their data, there is a hint of an increase in 
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Fig. 4. — Dependence of heliocentric mean velocity on metallicity (lower 
panel) and radial velocity dispersion (upper panel) for stars in our field. See 
legend in figure to identify the samples. The bins are chosen to contain the 
same number of stars, sorted by metal abundance; the horizontal error bars 
represent the range of metallicities in each bin. The vertical error bars are the 
eiTors on mean radial velocity and velocity dispersion, as appropriate. 

the velocity dispersion at l ower metal abundan ces, in agree- 
ment with our observations. IVieira et all (12007b find a flat dis- 
tribution of metallicity with velocity dispersion (from proper 
motion data) for stars with [Fe/H] > — 1, which is not in dis- 
agreement with our observations, where most of the increase 
in velocity dispersion takes place for lower metallicity stars. 

For stars with [Fe/H] > —1 the velocity dispersion is in 
good agreement with that measured for BRAVA M giants 
and does not depend strongly on metal a bundance, which is 
broad agreement with what measured by IVieira et al.l (l2007l) 
in Plaut's Window and th e two lower galactic latitude fields 
in iBabusiaux et all (120101) . However, at [Fe/H] < -1 there 
appears to be a (at face value) smooth transition to a dynam- 
ically hot component. Similarly, in the upper panel of Figure 
4 we see that the metal-rich component appears to have sig- 
nificant mean heliocentric velocity, while at [Fe/H] < — 1 one 
observes a smooth trend towards a relatively static velocity 
component. One caveat in this is that the metallicity errors 
are large, and we cannot rule our a bimodal distribution with 
the 'wings' of the errors simulating a smoother transition be- 
tween the two behaviors, although this would require a cor- 
relation between metallicity errors, measured radial velocity 
and velocity dispersion. 

The metal rich stars are best interpreted as part of the 
bar/bulge structure. Their kinematics show evidence of ro- 
tational support and bulk rotation and are consistent with data 
from the BRAVA survey in this region. The behavior of the 



more metal poor component, showing high velocity disper- 
sion, and low to zero velocity relative to the Sun may be ex- 
plained by a classical bulge or by inner halo stars. With a 
mean metal abundance of [Fe/H] —1.5 these stars appear 
to be best interpreted (at least provisionally) as an inner halo 
population, although bulges can of course be metal poor as 
well. T his i s consistent with t he ob servations bv lZoccali et alj 
(120081) and IBabusiaux et alj (120101) . albeit for a single sight- 
line. However, the BRAVA data show no classical bulge com- 
ponent fitting their dynamical model (IShen et al.ll2010i) . One 
possibility is that by selecting M giants and using the Calcium 
Triplet as their main radial velocity indicator, BRAVA may be 
biased against lower metallicity stars and therefore preferen- 
tially miss the high a component. 

The properties of galactic globular clusters present an inter- 
esting analogy with what is observed here: metal rich clusters, 
with mean [Fe/H] of ^ —0.7 are believed to be associated 
with the bulge and are supported at least in part by rotation, 
whereas inner halo clusters have mean [Fe/H] of ^ —1.6 their 
kinematics are dominated by random motions and at most 
very slow rotation. It is tempting to speculate that the two 
components we see in our data are analogous to the metal- 
poor and metal-rich globular cluste rs, wh ose properties they 
appear to share to some extent ('cf.. IBabusiaux et al. 2010 for 
a simi la r two-component mo del f or the bulge) . Ortolani et all 
( [19951) : IZoccah et alj (l2003l) and lClarkson et al. (2008) have 
argued, on the basis of isochrone fits to bulge globular clusters 
and field stars, that the bulge formed nearly coevally with the 
halo. Most globular cl usters in the inner halo formed within 
±1 Gyr of each other (iMarin-Franch et al. 2009) . If this ap- 
plies to bulge stars as well, it would imply a rapid star forma- 
tion process, at least for the inner regions ( < 20 kpc) of the 
Milky Way 

If this is the case, the smooth transition between the metal- 
rich and metal-poor subsystems, with a 'turn-over' point at 
[Fe/H] ^ — 1 may imply that the bulge and halo components 
are continuous and that there is no clear dichotomy between 
the two (modulo the large errors in metal abundance). This 
would be consistent with the BRAVA result that the bulge was 
formed (in a dynamical sense) from secular evolution at high 
redshift. As long as the stars also formed rapidly, the ob- 
served a-element enhancements are not in disagreement with 
this hypothesis. 
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